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The relation between the Seiberg-Witten prepotentials, Nekrasov functions and 
matrix models is discussed. The matrix models of Eguchi-Yang type are derived 
quasiclassically as describing the instantonic contribution to the deformed partition 
functions of supersymmetric gauge theories. The constructed explicitly exact so- 
lution for the case of conformal four- dimensional theory is studied in detail, and 
some aspects of its relation with the recently proposed logarithmic beta-ensembles 
are considered. We discuss also the "quantization" of this picture in terms of two- 
dimensional conformal theory with extended symmetry, and stress its difference 
from common picture of perturbative expansion a la matrix models. Instead, the 
representation for Nekrasov functions in terms of conformal blocks or Whittaker 
vector suggests some nontrivial relation with Teichmiiller spaces and quantum in- 
tegrable systems. 



1 Introduction 

The Seiberg-Witten (SW) prepotentials play an important role in studying properties of the 
gauge theories at strong coupling [1, 2] (see also [3] and references therein for recent discussion 
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of these issues). Simultaneously, they are beautiful objects of pure interest for the mathemati- 
cal physics, and have been recently re-derived as quasiclassical limit of the Nekrasov instanton 
partition functions [4, 5, 6, 7, 8, 11, 12]. The latter ones have the form of statistical models, 
where summing is taken over the sets of random partitions or Young diagrams [4, 13]. Originally 
arisen from the integrals over the instanton moduli spaces in the double-deformed supersym- 
metric gauge theories, they are related directly to the correlation functions or conformal blocks 
in two-dimensional conformal quantum field theories [14] (see also discussion of this relation 
e.g. in [15]). 

Quasiclassics of Nekrasov functions produce the SW geometry, so that the prepotential 
appears as a critical value of the free energy functional of the correspondent statistical model 
[7, 8]. This happens quite similar to the case of quasiclassics in matrix models, though some 
details - to be discussed below - in these two cases are still quite different. There have been 
already many suggestions (see recent attempts e.g. in [9, 10]) to treat the Nekrasov functions 
in a similar to the matrix models way beyond the quasiclassical approximation - at least in the 
sense of perturbative expansion, producing corrections to the prepotential. The relation seems 
however to be not very straightforward. 

We start the demonstration of this relation in the most old and transparent example. In fact 
it becomes more natural, if instead of full Nekrasov function one would take only its instantonic 
part, leaving aside the perturbative contribution. In the simplest case, one gets in such way 
the Eguchi-Yang (EY) one-matrix model [16], but in the limit of vanishing size of the matrix 
[7, 17]. The SW periods are introduced in this context in a completely different way, not being 
the fractions of condensed eigenvalues, like one used to have for the common matrix models. 

Partially this reflects the difference of how the string coupling - parameter of genus or 
quasiclassical expansion is introduced. For the matrix models it is related to the rank of 
gauge group, with sensible expansion only at its large values. In the picture of summation over 
partitions the string coupling is rather related to particular combination of the two deformation 
parameters, while the rank or matrix size tends to zero. Differently, the same distinction comes 
from the completely different role of the Virasoro constraints in these two theories. For the 
matrix models the Virasoro constraints are equivalent to the loop equations, giving directly the 
spectral curve and an iterative procedure of constructing the perturbative corrections, solving 
the loop equations order by order (this can be thought as one of the basic features of the B-type 
theory). In the case of Nekrasov functions, which rather belong to the A- type string models, 
the Virasoro constraints determine only the gravitational dressing of Nekrasov function by the 
descendants of unity operator, giving no restriction to the form of this function itself. 

We are going to discuss also the relation of SW geometry with the recently proposed logarith- 
mic matrix models [18]. For this purpose we study in detail the theories with the fundamental 
flavors and present, in particular, the explicit solution for the simplest of them - an Abelian the- 
ory with two flavors. It gives a hint, how the SW geometry can be rewritten in the form close to 
that of the matrix models, but this does not lead to the full identification of these two pictures. 
Instead, coming back to the pure gauge theories - where there is even no way for identification 
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of the spectral curves - we show how reformulation of the SW geometry suggests the way of its 
natural quantization from the angle of view of the relation with the two-dimensional conformal 
field theories, and this leads straightforwardly to a natural conjecture for the expression of 
Nekrasov function in terms of this quantized picture. 



2 SW prepotentials as limit of Nekrasov functions 

The instantonic calculus [4] in Af = 2 supersymmetric gauge theory gives rise [7, 8] to the 
extended SW prepotential as a critical value of the functional J 7 ~ — log Z, having an 

ei,2->0 6162 

integral representation 

T=\\ dxf(x)F uv (x) - \ f dx 1 dx 2 f"(x 1 )f"(x 2 )F(x 1 - x 2 ) (2.1) 

J ^ Jx 1 >x 2 

and giving the main contribution to the Nekrasov function Z(»|ei, e 2 ) at vanishing deformation 
parameters e± t2 — > 0, when extremized w.r.t. the second derivative of the profile function 
f"(x) = ^ of the Young diagram. In (2.1) the bare UV potential equals 

Fuvix) = Y,^TTi ( 2 - 2 ) 

fc>0 

and the kernel 

W = y (logx-0 (2.3) 

comes up from the (generalized) Plancherel measure in the sum over random partitions [6]. It 
is expected, that the main contribution at ei j2 — > is given by some large "limiting" Young 
diagram with the profile f(x), to be found by solving the variational problem for (2.1) upon 
normalization conditions imposed as constraints 

a* = \ I dx xf"{x), i = l,...,N (2.4) 

which can be in standard way taken into account by adding them to the functional with the 
Lagrange multipliers 

T -> T + a i (^i - \£ dx x f"( x )^j (2-5) 

(for the case of U(N) gauge theory one has to consider solutions with cuts {Ij}, % — 1 . . . , N). 

The whole setup for (2.1)-(2.5) is almost identical [19] to the standard quasiclassics of the 
matrix models (see e.g. [20, 21, 22, 23, 24, 25]), but with few crucial distinctions: 
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• The Coulomb gas kernel in (2.1) is replaces by a multivalued kernel (2.3). 



• The properties of the double derivative f"(x) of the shape function for the large extremal 
Young diagram are essentially different from those of the eigenvalue density in matrix 
models. 

The extremal equation for the (2.1) gives the system of N integral equations 

J2t k z k - dxf"(x)(z-x) (log \z-x\ - 1) = af, z eli, i = l,...,N (2.6) 

k>0 J 

on each segment of the support. Generally the solution can be expressed in terms of the Abelian 
integrals on the double cover 

y 2 = U(z-xt)(z-x-) (2.7) 

i=\ 

which is a hyperelliptic curve of genus g = N — 1. Define 

S(z) = F{j V {z) - J dxf"(x)(z - x)(\og(z - x )-\)-a D = 

n (2.8) 
= y^t k z k - 2N ■ z(logz - 1) + 2V cii • logz + . . . 

where the integral is taken over the whole support I = U^Ij, a D = Y^j=i a f \ an d consider 
its differential, or 



fc>0 

satisfying 



(2.9) 



<£>(x + iO) + $(x-iO) = 0, xelj, i = l,...,N (2.10) 
on each cut, and normalized in the following way 

$(4) = 0, 

Q(xJ±iO) = $(xf_ 1 ±iO) = ±2m(N-j + l), j = 2,...,N (2.11) 

®(xi) = ±2niN 

If all tk = for k ^ 1 and e* 1 = A 2N , the derivative <3> = ^ is an Abelian integral on the curve 
(2.7) with the asymptotic 

$ ^ T 2iVlog2±2A^logA + 0(^ 1 ) (2.12) 
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( Z (P±) — °°)) whose jumps are integer-valued due to (2.11), or § rf$ ~ 47riZ. It means that 
the hyperelliptic curve (2.7) can be seen also as an algebraic Riemann surface for the function 
w = exp (— $/2), satisfying quadratic equation 

AN ( w+ i) = Pn{z) = W z - w <) ( 2 - 13 ) 

^ ' i=i 

or 

y 2 = p N ( z f - 4A 27V (2.14) 
i.e. the branch points {xf} are roots of Pn(z) =F = 0, and 

y = A" («,-!) (2.15) 

Note, that in the simplest N = 1 case the form (2.13), i.e. 

z-u = A^iu + ^ (2.16) 

can be always achieved by change of the variables, i.e. it can be used as well for the switched 
on higher flows [7, 8]. However generally, for nonvanishing higher couplings in the UV potential 
(2.2) the profile function cannot be obtained as a jump of any algebraic function on the curve 
(2.7), unlike the case of resolvent and density in matrix model. 

The generating differential (2.9) is now 

dS= -2\ogwdz = -d(2z\ogw) + 2z— (2.17) 

w 

just the Legendre transform of the SW differential dSgw ~ z ^ on the curve (2.13), (2.14). It 
periods 

1 f dw 



<H = — <b z— (2.18) 
2m J A . w 

coincide with the SW integrals and the only nontrivial residues at infinity give 

res P+ (z^dS) = -res P _ (z^dS) = logA 2Ar 

N 

res P+ (dS) = -res P _ 



V 3 



(2.19) 



The differential (2.17) satisfies the condition 



r iri 5w , SP(z) , 

bdb ~ — dz = — dz = holomorphic (2 20) 

w V £f=i«i=° 



where the variation is taken at constant co-ordinate z and constant scale factor A. This provides 
the integrability of the gradient formulas 




(2.21) 



reducing it consistency to the symmetricity of the period matrix - a particular case of the 
Riemann bilinear relations. 



3 Eguchi-Yang matrix model 

Let us now substitute into (2.1) the shifted profile function [7, 17] 

f(x) = \x-a\ +g(x) 
f'(x) = sign(x — a) + g'{x) = sign(x — a) + p(x) (3.1) 
f"(x)=25(x-a)+g"(x) 

where the function g(x) and its derivative g'{x) = p(x) vanish at the ends of the cut I = (x-,x+): 

g(x±) = 0, g\x±) = p{x±) = (3.2) 

One obviously gets 

T = F uv {a) + \ [ dxg"(x)F uv (x) - [ dxg"{x)F{x - a)- 

1 3 3 (3.3) 



where 



~ I dxidx 2 g"{xi)g" (x 2 )F(xi - x 2 ) = F uv (a) - J" inst 

J X\>X2 

^inst = \ I dxp(x)W{x) - ]- I dxidx 2 p(x 1 )p(x 2 ) log(x 1 - x 2 ) 

W(x) = F[j V (x) - 2(x - a) (\og(x - a) - I) 
The variational problem for (3.4) is solved under additional constraint, following from (2.4) 

— J dx xg"{x) = J dxp(x) = (3.5) 

The dependence of J-i ns t upon the variable a can be easily obtained by redefinition of the 
re-variable x — > x + a and renaming p(x + a) — > p(x) 

^inst = 9 / dxp(x) (F' uv (x + a) - 2a;(loga; - 1)) - 

1 J (3.6) 
dx 1 dx 2 p(x 1 )p(x 2 ) log(xi - x 2 ) 



Xl>X2 
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so that 



and it means that 



= \ j dxp(x)F^ v (x + a) = J2 Hfe ^7 (3-7) 

k>l 



^inst(a; t\, t 2 , ■ ■ •) — J 7 o(tl, t 2 , • ■ ■; 



ti-H^ v (a) 



and 



^0 = -^inst I o=0 = \ j dx P( X ) ( F Uv( X ) - 2x(logX - 1)) - 

dx 1 dx 2 p(xi)p(x 2 ) log(xi - x 2 ) 



(3.8) 



(3.9) 



I Xl>X2 

is just the effective potential of the EY matrix model [16], constrained by vanishing of the total 
density of the eigenvalues (3.5). 

It means, that the dependence on zero Toda time is introduced here in quite uncommon for 
the matrix models way (3.7), (3.8). On the level of the Toda chain hierarchy, which governs 
the dynamics over the parameters of the potential [29], one deals here with two completely 
different solutions to the Toda equation 

namely 



exp ( ) (3.10) 



mamo 



= \t\ flogt - ^ + itof? + ••., t = hN (3.11) 

for the standard matrix model (see e.g. [25]), but 

^ey = \Ut\ + e' 1 + . . . , t = a (3.12) 

for the EY model. The essential difference comes in the zero-time dependence: for the standard 
matrix model free energy necessarily contains the function (2.3) (F(t ) in the r.h.s. of (3.11)), 
giving rise to the logarithm into the second to- derivative, "canceled" further by exponentiating 
in (3.10), while in the EY case the absence of logarithms of to — & requires necessarily the 
exponential dependence on t\. Moreover, the W-hosoxv masses, which are associated to the 
periods of dual to (2.8) SW differential, are not related with the filling fractions of the matrix 
model with the EY potential. 

Note also, since 

W(x) = F{j V {x) + 2a(loga - 1) + £ (3.13) 

n>0 V > 
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the effective potential of the EY matrix model satisfies the property 

J~inst — J~mst{tl', t 2 , ^3, • • •) 

f , (3.14) 

'' = i ' + tFi^r - k ^ 2 

and has a natural expansion at large values of a. 

For the non-Abelian case instead of (3.1) one has to make a substitution 

N 

f(x) = L(x; a) + g(x) = ^ \x - aj\ + g(x) 

3=1 

N N 

f'( x ) = $Z si & n ( x ~ a i) + = ^2 sign(x - Uj) + p(x) (3.15) 

N 

/"(*) = 2 + </'(*) 

3=1 

after which again the function g(x) and its derivative g'(x) = p(x) vanish at the ends of all cuts 
I, = (*-,*+), I = Uf =1 I,: 

g(xf) = 0, g'( x f) = p(xf) = 0, j = l,...,N (3.16) 

Now, instead of (3.3) one gets 

N 

F=J2 F uv(aj) - F ( a i ~ a j) + 

3=1 i+3 

N 



(3.17) 



+ i I dxg"(x)F uv {x)- I dxg"{x)J2F(x- aj )- 
J J j= i 

-- \ dx 1 dx 2 g"{x 1 )g"(x 2 )F(x 1 - x 2 ) = 

^ J Xl>X2 

N 

= ^2 F Uv( a j) — ^2 F ( ai ~ ~~ ^ inst = -^0 — ^inst 
3=1 i¥=3 

where F Q is just the sum of the classical and perturbative contributions to the prepotential 
(see Appendix A), while the instantonic part is again given by effective potential of a "matrix 
model" 

^inst = ~ / dxp(x)W(x) - - / dx 1 dx 2 p(x 1 )p(x 2 )\og(x 1 - x 2 ) 

^ J ^ Jx 1 >X2 

n (3.18) 
W(x) = F[j V (x) - 2 J> - a,) (\og(x - a,) - 1) 

3=1 



The variational problem for the functional (3.18) is again solved under constraint (3.5) at each 
component of the cut I, since 

J2 a ?( a i-l [ dxxf"(x) \ =-±J2 a f I dxxg"(x) = ±J2 a f \ dx P( x ) ( 3 - 19 ) 

We have found therefore, that in the non-Abelian case the instantonic partition function is 
described by the EY type matrix model with vanishing filling fractions, where the role of the 
SW periods is played by the impurities in the potential. 



4 Supersymmetric QCD and matrix models 

In order to analyze possible outcome of this similarity, let us now turn to the case of supersym- 
metric QCD, or the M = 2 supersymmetric Yang-Mills theory with extra fundamental flavors 
of matter [26, 27]. We shall concentrate mostly on four- dimensional conformal theory, i.e. with 
the number of flavors Nf = 2N and the vanishing J\f = 2 beta-function (3 = 2N — Nf — 0, 
coming back later to the case of pure supersymmetric gauge theories, after taking the limit of 
large masses of the fundamental multiplets. 



4.1 Abelian theory with two flavors: explicit solution 

In the case of Abelian theory with two flavors f — 1,2 the functional (2.1) is changed just by 
substitution Fuv{x) — > Fjjv(x) F( x ~ m /)> ( see Appendix A, where this is discussed from 
the point of view of perturbative prepotentials) 



2 

so that the function 



F=\J dxf'{x){F uv {x) + Y^F{x-m f ))- 
\ ( dx 1 dx 2 f"(x 1 )f"(x 2 )F(x 1 - x 2 ) + a D ( a - \ [ dx xf"[ 

^ Jxi>X 2 V J 



X 



S(z) = (z-m f )\og(z-m f ) + ..., f = 1,2 

z—trrif 



(4.1) 



S(z) = ^zW(z) = F{jv{z) + ^ F ' {Z " m/) " / F ' {Z ~ x ^"^ dx - ° D = 
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acquires extra singularities at 2— > mf, f — 1,2, while the EY z log z- singularity at z — > 00 is 
canceled due to vanishing beta-function (for two flavors). 

More transparently it is seen for 



*(z) = ^ = F^ v (z) + hg(z -m f )-J \og(z - x)f"(x)d. 
= khz"' 1 + £ log (l - ^f) - I log (l - X -) f"{x)dx 

fc>0 k>0 \ f 

$(z) = \og(z-m f ) + 0(l), f = l,2 



(4.3) 



and 



d * m , , , v- 1 r f"( x ) dx 



E— -( 

z — mf J 



, FZ v (z) + z 

dz z — mf J z — x 



= ^^~l)M fc - 2 - 2Q ^ fmf +o(^) (4-4) 

k>l 



dz 

d^{z) = + ..., / = 1,2 

z-^mf z — rrif 

The solution can be still constructed as an odd under y -b- — y differential on a cylinder 

y 2 = (z — x + )(z — X-) (4.5) 

d * = 1 ^ \ < 4 ' 6 > 

(z - mi){z - m 2 )y 

On small phase space, where only t\ 7^ 0, the asymptotic (4.4) requires numerator of (4.6) to 
be a linear function 

Hz)\ tk =t lSk>1 = ^M + ^o (4.7) 

where 

^1 = yi + V2 = ^2 m f ~ 2a 

f (4.8) 
V>o = -m 1 y 2 - m 2 yi 

and 

y f = y(m f ), f = 1,2 (4.9) 
The solution (4.6) can be conveniently described by the following anzatz 

d ^ _ 1 dz { 2v ~ g/ m f) z + 2m ^ - v J2f m f = (4 10) 

v^i - 4C 2 y (z-mi)(z-m 2 ) ~ 
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where we have introduced 



with 1 



W + ^ = ^=y Y> = (I - 4C 2 )y 2 = (z - v f - ±Q(z) 
Q(z) = C\z-m 1 )(z-m 2 ), y, = ™f^- , / = 1,2 



C= a/1 -4C 2 = -tanh- 

S 2 cosh |' v s 2 

^7==^ = - ^"5Z m /) coth| = 2a-^m / 



(4.11) 



(4.12) 



At e' 1 — > the last relation turns into v — a, and in this limit at large masses one can introduce 
finite scale C^m\m 2 = A 2 . 

In order to write explicitly the generation function (4.2) one needs to introduce first the 
(similar to (2.16)) uniformization of (4.5) via 

L ( 1 \ L ( 1 \ dm dz 

zu~\ = z - V, — I w \ =y, 



2 \ w J 2 \ zu J vo y 

= x + + x_ v - 2( 2 (m 1 + m 2 ) ^ . 

2 1-4C 2 
^ = x + _ x - = _ ^^ 2 a/C 2 (^i - rn 2 ) 2 + (v- m^iv - m 2 ) 



and define the auxiliary functions 



zu-zu f 
Xf-— 7, J -1,2 

ZDVJf — 1 

2 2(m f -\/)(^-\/)-L 2 ... 
X-?- W + l=0, = *_ m/ (4.14) 

Xf = Wf = y(m f -V + y f ), f = l,2 



The generation function (4.2) now reads 



S = ^2t k tt k + (2a - m 1 -m 2 )\ogzu - ^(z - m f ) log*/ - 2/logH 2 

k>0 f 

Q k = z k (w) + - z h (zu)- , k>0 (4.15) 
C = log(roiro 2 ) = ^ 



1 The sign of the root is chosen to fit with the weak-coupling regime, when e* 1 — > and tanh < 0. 
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where we have just used the basis of odd under the involution functions on the cylinder (4.5) 
with the only possible singularities at two infinities, e.g. such as 

tt = \ogzu, Sl! = y, ... (4.16) 

If only ti 0, 2 2 = e* 1 , and one gets from (4.15) 

S (*) lt fc =ti« fc ,i = ~ 2z lo S W + (2a-m 1 -m 2 )\ogzu+ m f lo S Xf (4.17) 

/=1,2 

where M^ 2 = X1X2 is defined in (4.11). 

Using these formulas it is easy to compute the resulting prepotential for (4.10), (4.12), which 
reads (for the only nonvanishing t±, and up to the linear terms ~ (mi + m 2 )a, which do not 
influence onto the second derivatives - coupling constants, and can be eliminated by adding the 
linear terms ~ (m x + m 2 )x to the potential in (4.1)) 

J 7 = \aH x - [a- m 1 )(a - m 2 ) log (l - e' 1 ) + J r pcrt (a; m) (4.18) 



and contains 



so that 



J r pert (a; m) = ^ F ( a ~ m f) 
f 



(4.19) 



d 2 T e' 1 \\f{a - mi) 

r{a) = ^ = log L ^ f> = log a N ' + W (4.20) 



where 



W = log — = log , (4.21) 

(l-e*i) 2 m/-*o & (l_ e *i) 2 

does not depend at vanishing masses on the vacuum expectation values (condensates) and gives 
rise to the non-perturbative renormalization of the coupling 

e T / 2 = - or 

e -ti/2 _ gti/2 ' 

^ e fc *i ( 4 - 22 ) 
r = t 1 -21og(l-e* 1 ) =tx + 2^ — 

fc>0 

which is a toy- model analog of the instanton renormalization of the coupling constant, given 
by the Zamolodchikov asymptotic formula [28]. 

The formulas for the EY matrix model (3.4) in conformal case remain almost intact, except 
for the potential 

x — mf) — 2F'(x — a) = 

_ f (4.23) 

= 2_.tkX k + / X x ~ m f) Q-°s( x ~ m f) ~ 1) ~~ 2(x — a) (log(x — a) — 1) 

fc>0 / 
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For a = and m/ = 0, f — 1,2 the potential (4.23) turns back into the potential a standard 
1-matrix model. 



4.2 Different parameterizations of the curves 

We have already seen in previous section, that to construct the exact solution explicitly one 
rather needs the parametrization (4.13), than more common for supersymmetric QCD form of 
the curve (4.11). Let us point out now, that the curve (4.11), can be also presented in the form 

qi(z)w 2 - (z - v)w + q 2 (z) = 

Qf(z) = C(z-m f ), f = l,2 ( 424 ) 



2 ^2 TJ T 2 12 

w = —W = —X1X2 
Qi Qi 



Introducing another new variable 



x _ _ z /<M^ 



w W y q 2 (z) 
one can rewrite (4.24) as 

(mxw 2 - vw + (m 2 miw 2 - v (S + |) w + m 2 



(4.25; 



x = 



(4.26) 



m 2 a — m 2 a — m\ 
w w — S w — i 



where we have used (4.12), (remind also that on small phase space S 2 = e* 1 is just the expo- 
nentiated UV bare coupling). The SW differential becomes 

dW dw z ( dqi dq 2 \ , m\ dz m 2 dz 

dSsw = z-— = z h = xdw + — = 

W w 2 \ qi q 2 J 2 z - m 1 2 z - m 2 f4 27) 

in , mi dz m 2 dz 

= dS G + 

2 z — mi 2 z — m 2 

The first term in the r.h.s. dSc = xdw, and exactly in such form has been introduced in [31], 
as appearing naturally in the context of D-brane considerations, it is normalized due to (4.26) 

as 

ies w=0 xdw = m 2 , res w=OQ xdw = —mi 
res w= zxdw = a — m 2 , res w= i xdw = —(a — mi) ^ ^ 
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This form of the curve is especially natural from the point of view of comparison with the 
conformal representation [5] in terms of a theory of two-dimensional free scalar field 



Z(x;a, mi ,m 2 ) = Z pcTt ■ ( e ^i^(oo) e i(a-m 1 )^(l) a .L 0e -i(a-m 2 )^(l) e -im 2 ^(0)^ = 

2 2 

= Zp Crt • x^ ^™^ ^ e <"*i0(oo) e i(o-mi)^(l) e -i(o-m2)»(») e -im20(O)^ = (4.29) 

= Z pert -x a2/2 (l-x)- {a - mi)i " 2) 

of the U(l) partition function with Nf = 2. One can absorb the perturbative contribution 
inside the correlator by redefinition of the scalar product in two-dimensional theory, see the 
details below. 

From (4.26) one also easily finds, that 



d dx (5 - h) dw 
— — don = —aw = — ; ^-r- = du 

da G da ( w _ S )( w _i) ( 4 . 30 ) 

res^s du> — 1, res„, = ^ du = — 1 



is just a degenerate analog of the holomorphic differential. 



4.3 Supersymmetric QCD, conformal theory and logarithmic po- 
tentials 

In the case of nonabelian supersymmetric QCD (the TV-cut solution) the formula for the EY 
matrix model potential (with fundamental matter) obviously changes for 

N f N 

W(x) = F' uv {x) + ][> " m f ) (log(x -m f )-l)-2 J> - a,) (\og(x - a,) - 1) (4.31) 

f=l 3=1 

It means, that in conformal case Nf = 2N the derivative of this potential does not have any 
longer a logarithmic singularity at infinity 

W'{x) = F^ v {x) + - m f)- 2 Yl lo ^ x ~ a ^ = 

(4.32) 



F» vW + ^l„ g (l-^)-2i:iog(l-|) 



but still has logarithmic singularities at x — aj and x — rrif. 

Generally, for the U(N) theory with Nf < 2N flavors on small phase space one usually 
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starts [26, 27] with the analog of the representation (4.11) 

N-2 N f (4.33) 

p(z) = z n ~Y1 Q( z ) = A2N ~ Nf U( z - m /) 

k=0 f=l 

with the generating differential 

^H p ' {z) - PJ wf) (4 ' 34) 

In the conformal case Nf = 2N instead of the scale parameter A one gets the dimensionless 
constant, and equation (4.33) is changed for 



W + ± = -^=, Y 2 = (1 - 4C V = P(zf - 4Q(z) 

2N 

Q{z) = C\[{z-m f ) 



(4.35) 



For the vanishing masses the differential (4.34) becomes holomorphic 

dSsw = N % (zP'(z) - ^fp(z)) = 

Q(z)=cz N f y v 2 j 



Q(z)=C,z 

dz I / N f \ „ ( N 



nl *-?)^E(?-')^)= («6) 



Nt=2N Y 
1 k=0 



For example, in the N = 2, Nf = 4 case the differential (4.36) is proportional to the only 
holomorphic differential on torus (see e.g. [28]) 

,_, dz u dz d^ 

dS sw ~ u— = (4.37) 

y y/i - 4C 2 y v 



with 



v 2 = n (e-g/)=^-i)(e-A) 

/=0,1,A 

2 A 1 1 (4-38) 

c - 
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and since S = e T0//2 this is identical to the formula (4.12) from the N — 1 case. 

However, for the nonvanishing masses it is more convenient to use the analog of the repre- 
sentation (4.24) 

q(z)w 2 - P(z)w + q{z) = (4.39) 

where P(z) is the same polynomial of power N as in (4.33), while for Nf = 2N one can naturally 
introduce 

N N 

Q(z) = C ~ m f ), q(z) = ( JJ(z - m f ) (4.40) 
/=i /=i 

for some natural splitting of the Nf = 2N matter multiplets into N "fundamental" and N 
"anti-fundamental" - with the masses m and m correspondingly. Again, like in (4.24), we have 
introduced here 

w * = M W 2 (4 41) 

q(z) 

For iV = 2 in (4.40), equation (4.39) becomes as well a quadratic equation in z- variable, and 
can be therefore written [321 as 



C 2 (w)z 2 -C^z + Coiw) = (4.42) 

with 

C 2 (w) = (w 2 - w + C ^ = i (e + (w - E) (^w - (4.43) 
Equation (4.42) can be further re-written as 

2 



^I'-sU =ss-S !s * v (444) 



2 _ / z Cl \ _ Cl Co ,2„,.2 



2 / <-v2 

or 

d(tt;) 2 C ( W 

4-u7 2 C 2 (u>) 2 w 2 C 2 (w) 
At vanishing masses the last equation turns into 

2 _ Cq(w) m m 1 

x = ~w*c 2 (w) = ~cw( W -e) ( w -i) = ~c^m-m-\) (4,46) 

in the notations of the elliptic curve (4.38). 

Generally, instead of the holomorphic differential (4.37) one can write 
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for some polynomial 



M0 = 5>e ^ n C -A) 

j=0 



(4.48) 



reproducing the holomorphic differential (4.37) in the limit of vanishing masses. Due to relation 
between the coefficients of and the residues 



res ?=0 dS S w = y/pifi) = y/Po = Am 

res f=oc dS S w = y/pl = -m^ 
res f= i dS S w = VpU) = (1 - A)mi 
res ?=A dSW = Vp(A) = A(A - l)m A 



one gets 



£moduliP(0 = ^moduliC ' £(£ ~ !)(£ ~ A) 

and it means that for the variation of the SW differential one can write 



^moduli dSs W 



•'moduli 



p(0 



VROe(e-i)(e-A) 

so that the r.h.s. in this equality is obviously holomorphic. 
One can also present the 2-differential T(d£) 2 = (dS S w) 2 a s 



moduli 



c 



di 



r 



dS, 



SW 



/=0,1,A 



m. 



+ 



d£ 
Cf 



p(0 



e(e-l) 2 (e-A)^ 



(£-<7/) 2 



- E 



+ 



c(0 



v/=0. 



e(e-i)(e-A) 



where 



and 



^ C / = 0, Ci + AC A = - ]T 



/=0,1,A 



c(0 = 



\/=0,l,A / 



+ [C — 2rriQmi\ — 2m m A ] 



(4.49) 



(4.50) 



(4.51) 



(4.52) 



(4.53) 



(4.54) 



is a linear function satisfying 5 mo duiiC(0 = A5 mo d u iiCo. Formula (4.52) may be assigned with 
two possible interpretations: 

• It has a sense of average of the stress-energy tensor in some two-dimensional conformal 
field theory 

\11a=o,i,a,oo VaWa)) 
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with the insertions of four primary operators of dimensions m 2 A , A = 0, 1, Aoo (m^ can be 
determined from the second equation in (4.53)) and Cf being the corresponding accessor 
parameters (see e.g. [30] and references therein); 

• The r.h.s. (4.52) has an obvious form of the matrix model curve, if one takes formally in 
the matrix model the following logarithmic potential 

V(0 = J2 m f lo ^ _ 9/) (4.56) 

/=0,1,A 

and the SW periods can be then identified with the filling fractions [18] . 

We do not find the last observation to be very useful, despite of existing vast list of literature 
on the subject, and in the rest of the paper we would discuss rather the first one, digressing it 
back to the case of the pure supersymmetric gauge theory. 



5 Pure gauge theories and Whittaker vectors 

If one takes in (4.24) the decoupling matter limit ( — > and m 12 — > oo so that A 2 = ( 2 mim 2 = 
fixed 2 , one gets instead of (4.26) the equation 

x = — -\ h A = ' \ /' ' = (J(w)) 5.1 

w 2 w (*|^) 

which is an average of the [/(l)-current J(w) = Ylnei over the state, satisfying 

J n |tf) = 0, n > 1 
Ji|*)=A|*), J |*)=a|*) 

Since [J n , J m ] = n5 n+mt o, equations (5.3) can be immediately solved explicitly by 

I*) = A^e* 7 - 1 ^) 
J n \a) = 0, n > 0, Jo|«) = o|a) 

i.e. in terms of the coherent state in the (charged) Fock module, so that 

(tf |#) = (ale^A^e- 7 - 1 ^) = A a2 e A2 = exp (|a 2 ^ + e* 1 ) (5.4) 

A 2 =e'i 

The generating differential can be treated as just as an average of the ?7(l)-current 

(J(w)) = dSsw — %dw (5.5) 



(5.2) 



(5.3) 



2 Literally, one needs to perform this limit in symmetric way, putting A = — (mi — —(iri2- 
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and its more conventional form (2.16) is restored by substitution x — ^. 

For the pure U(2) supersymmetric gauge theory instead of (5.1) one gets an equation [31] 



, A 2 u A 2 (V\T(w)\y) 



x 2 = — + — + — = v / = (T(w)) (5.6) 

w 3 w 2 w (*|^) 

which has an obvious sense of averaging of the stress-tensor T(w) = Xlnez ~~^T2 over the Whit- 
taker [12, 33] state \^) G "Ha,c i n the Verma module of the Virasoro algebra with the highest 
weight A and the central charge c, satisfying L n \ty) = if n > 1 and 

Li|tf)=A 2 |tf) (5.7) 

An important consequence of (5.6), proven in [34], is that the Nekrasov function of the corre- 
sponding theory is given by the scalar product, or the matrix element 



where 

after identification 



Z{a-e u e 2 ) = (tf|tf> = (^[A 4 ^^) (5.8) 
l*i>= I*>Ia=i (5-9) 



+ T^— , c=l + — , e = ei + e 2 (5.10) 



eie 2 4eie 2 ' e x e 2 



Moreover 



(5.11) 



|*>| A=0 = |fi)G^ A , c 

L n \n) = o,n > o, L \n) = A\n) 

i.e. \fl) is proportional to the highest-weight vector, but with a non-standard normalization 



~ r 2 (a|ei, e 2 )r 2 (-a|ei, e 2 ) (5.12) 

to the product of the Barnes double-gamma functions [35] (some details and other references 
are collected in Appendix). 

The nontrivial part of (5.6) contains the only equality 

being a sort of renormalization group equation. Again, by x = ^ the curve (5.6) turns into 
that with N = 2 from the family (2.13), and the generating differential is now [31] 

^{T(w))dw = dS sw = xdw = z^- (5.14) 
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Generally for the pure U (N) supersymmetric gauge theory one should write the average of a 
formal differential operator 

{V\V N \V) = 

V N = D N - T(w)D N - 2 - W^(w)D N - 3 - ... - W^H (5 ' 15) 

with the coefficients acting in the highest-weight module "H aiC of the W^-algebra, and is 
now the Whittaker vector 



W[ N) \V) = A N \^) 
Wi N) \V) = 0, n > 1, Wi K) \*) = 0, n>0, K<N 



(5.16) 



in the module with the highest weight \Q) = |^)| A=0 



N ( a K \ 

W^IQ) = s K (a; ei, e 2 )\Q) = V \ -L- + . . . K = 1, . . . , N (5.17) 

j=i V 1 " 2 J 

where also everywhere 

D\*)=x\V) (5.18) 

This set of equations (again by substitution x — ^) gives rise to the SW curves (2.13) endowed 
with generating differential dSsw = 2^u> = z— , with a natural conjecture 

Z= {^\k 2NL °\^) (5.19) 

provided by normalization 

(n\n)= (*|*)| A=0 ~ l[r 2 (a i -a j \6 U 6 2 )~z peit (5>20) 

written again in terms of the Barnes double-gamma function. 

In the N — 1 and N = 2 cases equations (5.16) (i.e. (V\D - J(w)\V) = for W 1 = C/(l) 
and (fy\D 2 — T(w)\^) = for W 2 = Vir can be explicitly solved w.r.t. x = (D) variable, and 
we come back to (5.1) and (5.6) correspondingly. To switch on higher times one just has to 
generalize the Hamiltonian in (5.19) 

Z(t u ...,t N )~ exp ^2t K W^j |*> (5.21) 

by introducing higher W- flows. Arising of the Gelfand-Dikij-type operator (5.15) in this context 
together with the W-flows in (5.21) suggests a nontrivial connection with the so-called W- 
gravity [37] and geometry of generalized Teichmuller spaces. 
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6 Conclusion 



We have discussed in these notes some relations between the formulation of supersymmetric 
gauge theories, coming from the instanton partition functions, and the matrix models. As we 
pointed out - there are several parallels of this kind, but neither of them seems to be very 
essential by itself. 

Let us finish instead with the following important remark. The quasiclassical picture of the 
matrix model is indeed very similar to the SW theory, and they are both described in terms 
of a quasiclassical integrable system, so that the prepotential is in fact a restricted Krichever 
tau- function [38]. The simplest part of integrable dynamics in this case is parameterized by the 
parameters of potential (2.2), and reproduces in the simplest case the well-known hierarchy of 
the dispersionless Toda equation (3.10) (for the higher-genera analogs of dispersionless equa- 
tions see [39]). The dependence on smooth periods (2.18) is far more transcendental and, on 
quasiclassical level, the knowledge is basically exhausted by the gradient formulas (2.21) and 
their consequences, like residue formulas, the WDVV equations, etc (see e.g. [38, 23, 25]). 

The formulas of sect. 5 demonstrate, that the nontrivial integrable dynamics over quasiclas- 
sical period variables has even far less trivial dispersive analogs. While the naive "quantization" 
of the dynamics over the times (2.2) leads (in absence of the "smooth variables" (2.18)) just to 
the string solution of the full hierarchy of the Toda or KP type with the tau-functions presented 
by the matrix elements in the two-dimensional theory of free fermions, the dispersive analogs 
of the full Krichever tau-function are matrix elements in non-trivial two-dimensional conformal 
theory, generally with extended symmetry. They also seem to be directly related with the 
quantization of Teichmuller spaces [40] (higher Teichmiiller spaces for the case of ^-gravity), 
Liouville theory [30, 41], and quantum-mechanical integrable dynamics in the systems of Toda 
type (see e.g. [42, 43, 44]). All these relations clearly deserve further investigation, and we are 
going to return to them elsewhere. 
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Appendix 



A Perturbative prepotentials 

The instantonic expansion J 7 = Ylk>o m ^ ne n on-Abelian theory starts with the perturbative 
contribution (which does not take into account the topologically nontrivial configurations of 
the gauge fields) 

N 

defined entirely in terms of the ultraviolet or classical (2.2) and perturbative (2.3) prepotentials. 
It is totally characterized by degenerate differential of (2.9) 

d% = F»' v (z)dz - = F^ v {z)dz - 2 jr JL. (A.2) 

and the coefficients of the polynomial Pn(z) (A. 2) (the same as in (2.13)) coincide with the 
perturbative values of the SW periods 

a, = -±Tes Vi zd$ = Vi, i — 1, . . . , N (A. 3) 

The perturbative generating differential is dSo = $odz, with 

N 

% = F^ v (z)-2j2^g(z-v,) (A.4) 
j'=i 



and satisfies 



^ = 2^, J=l JV 

Odj Z — Vj 

ddS fc _! 



(A.5) 



= kz K - L dz, k > 



what gives rise to 



Equations (2.21) 



N 

S (z) = F' uv {z) - 2 J> - Wi ) (log(z - - 1) (A.6) 



«? = £r = 5 oK0 (A.7) 



^0 



(9a 



completely determine (A.l), since on the perturbative stage one makes no difference between 



Vj and a,j. 
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For the theory with Nf fundamental multiplets instead of the formula (A.l) one has 

N N N f jv 

J^O = + ^pert = Yl F Uv(aj) - F ^ ~ ^ + Yl Yl F ( a i + ( A - 8 ) 

i=i i^j /=i i=i 

which can be obtained from (A.l) just by formal modification of the UV prepotential via 

N f 

F uv (x)^F uv (x)+J2 F ( x + m f) ( A -9) 

/=i 

which can be further used, after its substitution to the functional (2.1) to compute the full 
partition function for the theory with matter. 

The "quantum" - or, better, double- deformed - version of the perturbative prepotentials 
(A.l), (A. 8) can be written in terms of the Barnes double-gamma functions 

r2(x|ei ' e2) ~ n x+ J 1+me2 (a.io) 

n,m>0 



where the infinite product can be understand, say, via the zeta-regularization (see e.g. [36]) 

(A.H) 



logr 2 (x|ei, e 2 ) = ^ ^2 (x + ne 1 +me 2 y 

n,m>0 



s=0 



analytically continued to s = 0. The relation to (A.l), (A. 8) is established via the asymptotic 

logr 2 (x|ei, e 2 ) = — F(x) + less singular = 



log x - - + 



ei e 2 2 V 4 
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